Non-mean-field Critical Exponent in a Mean-field Model : Dynamics versus Statistical Mechanics 
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The critical exponent of susceptibility is expected as 1 in systems belonging to the mean-field universality 
class. On the other hand, the mean-field nature permits to introduce the Vlasov equation in the place of the 
Liouville equation, and the linear response theory based on the Vlasov equation gives the strange exponent 1 /4 
in the low-temperature side of the Hamiltonian mean-field model. We clarify that this strange exponent is due 
to existence of Casimir invariants which traps the system in a quasistationary state for a time scale diverging 
with the system size. The theoretical prediction is numerically confirmed by A'-body simulations and numerical 
Vlasov solutions. 
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Is the universality of critical exponents extended from the 
context of statistical mechanics to of dynamics? For a ferro- 
magnetic model, we can define the isothermal susceptibility 
X^ and the static isolated susceptibility ;i^'. The former comes 
from the statistical mechanics of isothermal processes, and the 
latter comes from dynamical arguments and is defined by the 
linear response theory 11 \. These two susceptibilities satisfy 
the inequality x^ ^ Af ' 121 Isft which takes care of the exis- 
tence of invariants in the considering processes 141 15|]. The 
above inequality implies a similar inequality between the crit- 
ical exponents y^, since susceptibility is scaled as;^^ ~ kr''* 
around the critical temperature T^. Here t = (J - Tc)/Tc 
is the reduced temperature, T is temperature, and 7+ and y_ 
are defined in the high- and the low-temperature sides respec- 
tively. Breaking of the equality may happen between the value 
of the susceptibilities, that is x^ > X^^ but this breaking does 
not imply the breaking of equality between the corresponding 
critical exponents. 

To tackle this problem, we concentrate on systems in the 
mean-field universality class, whose critical exponents are 
y± - I. One of such systems is the Hamiltonian mean-field 
(HMF) model la-tSI, in which the mean-field approximation 
is exact in the limit of large number of particles. HMF is an 
important toy model in the studies of long-range interacting 
systems and, moreover, it is used in the theoretical description 
of collective atomic recoil laser experiments H . The potential 
term of the HMF model is the same of the globally-coupled 
XY model, and hence the HMF model is a ferromagnetic 
model belonging to the mean-filed universaUty class. Mag- 
netization M in the dynamical context, namely time average, 
agrees with one in the context of statistical mechanics when 
initial conditions for A^ particles are sampled from the canon- 
ical equilibrium distribution and the external magnetic field 
is not introduced JlOll . In other words, the initial question is 
positively answered in the case of the magnetization, which 
scales as M ~ {-rf with/? = 1/2. This agreement allows to 
consider the initial question for y-t, since, to derive the expo- 
nents, we have to prepare canonical equilibrium states before 
adding the external magnetic field. 

For the HMF model, we can compute x^ explicitly and its 



critical exponents are y^ = 1 11 1111 . In contrast, the computa- 
tion of ;^'' is not an easy task since we have to solve the A^-body 
dynamics. Nevertheless, the mean-field nature of the HMF 
model permits to compute a dynamical susceptibility, since it 
naturally introduces the Vlasov equation in the place of the 
Liouville equation. Thanks to integrability of the one-particle 
Hamiltonian associated with the Vlasov equation, the linear 
response theory based on the Vlasov equation gives the exact 
susceptibility in the mean-field limit. We call it the Vlasov 
susceptibility and denote with;^'^. 

The critical exponent of x"^ is tT = 1 in the high- 
temperature homogeneous side [!l2|, but seems smaller than 
1 in the low-temperature inhomogeneous side lll3ll . We will 
show that y^ - 1/4 due to existence of invariants, called 
Casimirs. We clarify that existence of the Casimirs traps 
the system in a quasistationary state (QSS) and induces the 
strange exponent. We also show that the life time of these 
QSSs is 0(N), and hence, the accessible critical exponent of 
the susceptibility is not 1 but 1 /4 if the number of particles is 
large enough. 

We remark on studies of the critical exponents in the Vlasov 
dynamical systems. Based on the theory constructing an un- 
stable manifold in the Vlasov-Poisson equation ll411 . Ivanov 
et al. lll5n have stated that scaling laws are different from ones 
predicted by the ordinary mean-field theory. However, they 
start from an unstable spatially homogeneous Maxwell dis- 
tribution, and have said nothing about the critical exponents 
starting from the canonical equilibrium states. We stress that 
a special choice of initial states may give non-mean-field crit- 
ical exponents, but our analysis is referred to canonical equi- 
librium states. 

The HMF model consists of A^ particles confined in the unit 
circle on the (x, y) plane, and is expressed by the Hamiltonian 



H^ 



N \ 2 N ^ 

i=l L ^ 7=1 



- cosiqi - qj) 

2A^^ 



- h@{t) cosiqi - (p) 



(1) 



where h and 4> are the modulus and the phase of the exter- 
nal magnetic vector {h cos 0, h sin (p) respectively, and 0(f) 



is the Heaviside step function. The magnetization vector 
«M,)w, (Myh) is defined by (M,)^ + /(M,)^ = Z^ e'^'IN, 
where (■)^ represents the average over A^ particles. This sys- 
tem has a rotational symmetry with h - 0, and we may set 
{My)n - without loss of generality. We assume (p - Q 
to concentrate our interests on susceptibility along the direc- 
tion of the spontaneous magnetization, defined x. The mag- 
netization observable is M = cosq, and its average over 
the one-particle distribution f{q, p, t) is denoted by (M) - 
ff Mfdqdp. The corresponding efi'ective one-particle Hamil- 
tonian is written as 



On the other hand, the Vlasov susceptibility x^ is obtained 
by introducing the angle-action variables (6, J) in the Hamilto- 
nian system ^ with /i = 111611 . Thanks to the Vlasov linear 
response theory [[131, the Vlasov susceptibility is given in a 
similar form 



X''^(l-D^)/D\ 



and 



D' 



{M\ - <C2)o 



(10) 



(11) 



H[f](q, p,t)-p 12- <M) cos q - h&(t) cos q, (2) where C is the average of cos q over the angle 0, 



where we denote (M) - ff Mfdqdp the average over the one- 
particle distribution function f{q,p,t). This distribution is a 
solution to the Vlasov equation 



d,f + d„Hlf]dgf - d,H[f]d„f = 0, 



(3) 



associated with the effective Hamiltonian (|2]). 

To clarify the discrepancy between the isothermal suscep- 
tibility and the Vlasov susceptibility, we first compute the 
isothermal susceptibility for the canonical equilibrium state 
with temperature T. The mean-field approximation, which is 
exact in the limit A^ ^ oo in the HMF model, allows to com- 
pute magnetization by the one-particle distribution function 



//. 



^-[p-/2-((M)i,+h)cosq]/T 



rg-[/>V2-«M)„+/,)cos,]/r^^^^ 



(4) 

where (M)/, is the average of M with respect to //, . We put 
h - Q for f < 0, and h > Q for f > 0. The equilibrium 
magnetization (M)/, must satisfy the self-consistent equation 

{M)h = h (({Mh + h)IT)lh i{{M)h + h)IT) , (5) 

where /„ is the modified Bessel function of the first kind. 

From the definition of the isothermal susceptibility x^ ^ the 
magnetization (M)/, is written as 



{M)h = <M)o + hx'' + 0{h\ 
and hence //, is expanded as 



fh 



1 + h- 



-{M - (M>o) 



/o + Oih'). 



(6) 



(7) 



Multiplying //, by M and integrating over q and p, and com- 
paring the obtained expression with (|6]l, the isothermal sus- 
ceptibility is written in the form 



where 



X^ ^(l-D^)/D^ 



^^^^_ {M\,-{MY, 



(8) 



(9) 



C{J) 



- 1. c 

27r Jo 



cos q(9, J)dO, 



(12) 



and (C^)o = Jf C^fodOdJ. We compare the scaling of the 
Vlasov susceptibility ( fTOl ) with the one of the isothermal sus- 
ceptibility dH), both in the high- and the low-temperature sides. 
In the high-temperature side we have {M)q - and C = 0, 
since the angle variable is q itself, and (M^)o = 1/2 = T^. We 
therefore have D^ ^ D'^ ^ (T - T^)/T and 



X =X 



(13) 



which implies that the critical exponents are the same and of 
the mean-field universality class: 7+ - yX - ^- The theoret- 
ical prediction for x'^ is confirmed by direct numerical simu- 
lations of the Vlasov equation for small enough external field 
h as shown in Fig.[il The discrepancy between the canonical 
prediction and the numerics in Fig. [T]comes from nonlinear ef- 
fects. No QSS features can be seen in a small h region where 
the linear response theory is valid. 

In the low-temperature side, we separately analyze the com- 
mon term 1 - (M^)o/r and the individual terms between 
D^ and Z)^. The common part approaches to zero linearly 
again in the limit t /■ 0. Indeed, the mean square of M is 
{M-}o = 1 - r. Thus, we find 1 - {M\)/T ~ 2t. Moreover, 
the individual term of the isothermal function D^ ^ shows 
the scaling (M)^ ~ (-t)^^ with/3 =1/2 of the mean-field uni- 
versality class. The isothermal susceptibility, therefore, shows 
the classical exponent yl = 1 as expected by the statistical 
mechanics of mean-field models. On the Vlasov counterpart, 
by using the new variable k - -JjHlfo] + (M)n)/(2<M)n), the 
individual term (C)q in D^ (fTTl i can be expressed as 11311 



(C')o = 16 
-hl6 



Jo 



2E(k) 
Tik) 



kK(k)Mk)dk 



Xa 



2k^E{\lk) 

m/k) 



-2k^ + l 



K(l/k)Mk)dk, 
(14) 



where foik) = Aoexp[(M>o(l - 2k^)/T] and Aq is the nor- 
malization factor We note that the two integrals of (fT4l i are 
positive and finite. In the limit of ^ y 1, the divergence of 
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FIG. 1: (color online) Short time average of magnetization, M, di- 
vided by h in the high-temperature side. T = 0.51 (red circle), 0.55 
(green triangle) and 0.60 (blue diamond). Averages are taken in the 
period 100 </ < 3000. Temporal evolutions are obtained by the 
semi-Lagrangian method lUTIl , and the grid sizes are 256 x 5 12 (filled 
symbols) and 128 x 256 (empty symbols) in (q, p) e {-n, n] x [-5, 5]. 
Time slice is 0.05. Horizontal lines mark theoretical predictions of 
susceptibility l ll3b for T = 0.51, 0.55 and 0.60 from top to bottom. 
The broken lines are the predictions given by the canonical ensemble 
for the same set of temperature. 
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FIG. 2: (color online) The Vlasov susceptibility ;^''^ (red solid lower) 
and the isothermal susceptibility x^ (green broken upper) as func- 
tions of -T in log-log plot. The adiabatic susceptibility x^ (blue 
dashed middle) is also reported, which is computed explicitly by 
using the exact solution in the microcanonical statistics |lg| and by 
taking the invariance of the entropy during the quasistatic adiabatic 
process into account. The black straight lines are guides of eyes for 
the slopes -1/4 and -1. Purple square points are computed by time 
averages over the period < f < 500 in /V-body simulations with 
N = 10* and time slice 0.1. For each T, 10 points are plotted corre- 
sponding to 10 realizations. 



K(k) is written as log(l - k') and is integrable. In a large 
k region of the second integral, damping of fo{k) becomes 
slower as r y Tc since (M)() \ 0. However, the factor 
of the brackets gives 0{k^'^) damping, then the second inte- 
gral converges. Consequently, (C^)o approaches to zero in 
the limit r /■ as {C\ ~ ^/{M}^ ~ (-rf^, where we 
used the agreement of {M)o between the contexts of statistical 
mechanics and of dynamics llioll . The function D^ is hence 
written in the form D^ = ci\t\ + C2|t|^''^, and the second term 
dominates around r < 0. Accordingly, the susceptibility x^ 
scales as x^ ~ (—tY^^'^ and the Vlasov critical exponent is 
y^ - y6/2 - 1 1 A. This exponent is numerically confirmed by 
computing the two integrals in (fl4l i and by A^-body simula- 
tions as shown in Fig. |2] 

The strange divergence of the Vlasov susceptibility comes 
from the term (C^)o, and we will show that this term comes 
from the existence of Casimir invariants. A Casimir is written 
by using a smooth function ^ as J s(f(q,p, t))dqdp and is an 
integral of the Vlasov equation even if the Hamiltonian //[/] 
depends on time. The distribution function with the external 
field must keep the Casimir invariants, thus it is no more //, 
but /c. Writing /c - fa + g we have the restriction for g so 
that every Casimirs are invariant up to the linear order We 
remark that /o depends on {q, p) solely through //[/o], and ac- 
cordingly, is a function of the action / only. Using the Fourier 
transform of g with respect to the angle 9, the above restriction 
is rewritten as 



0: 



// 



[s{fo+g)-s(fQ)]dqdp- 



/'■ 



(MJ))go(J)dJ (15) 



where go{J) is the Fourier zero-mode of g. This constraint 
must hold for any smooth functions s, and therefore, we con- 
elude that goiJ) - B181I . We remark that the normalization 



condition for /c automatically holds accordingly. 

The above observation constructs g = /c - /o by subtracting 
the Fourier zero-mode from /;, - /q. For this purpose, we ex- 
pand/;, as O but replace ;^'^ with;^'^. The Fourier zero-mode 
of //, - /o is expressed by 



1 r" 

^J ifh-fo)dO: 



h^^(C(J)-{M}o)fo + 0{h'-). (16) 



Consequently, the constructing distribution function /, is 

/c = /o + h^^(M - C(J))fo + O(h^). (17) 

The distribution function /^ gives magnetization (M)^ = 
Jf Mfcdqdp in the form 

(M), = (M>„ + h^^-4^{{M\ - {C\) + 0(h\ (18) 



This equation and the expansion {M)^ - {M)o + hx^ + 0{h^) 
gives the Vlasov susceptibility x^ ( flOl l and the function D^ 

(El). 

Following the above discussions, we propose a scenario of 



relaxation as follows II19LI20I] . For f > 0, the system is trapped 
in a QSS to keep the Casimir invariants and this gives the 
strange ciitical exponent y[ - \ /4. However, the Vlasov dy- 
namics is not the true dynamics for finite systems, thus the 
Casimirs are not exactly invariants but evolve in a time scale 
which diverges with A^. Consequently, the system goes to the 
thermal equilibrium recovering the classical exponent after a 
long time. 

The above scenario is examined by direct A^-body simula- 
tions in Fig. |3] For f < 0, the system is in thermal equilibrium 
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FIG. 3: (color online) Af-body simulations of the HMF model with 
external field. N = 10^ (100), 10-* (10) and 10' (1), where the in- 
side of braces are the number of realizations over which the orbits 
are averaged. The system is in thermal equilibrium with T = 0.499 
in f < 0, and the external field turns on with h = 0.01 at f = 0. The 
direction of magnetization vector (Af j, M,.) is reset to x-direction at 
f = 0. In each panel, three horizontal lines represent the thermal equi- 
librium level without external field (lower), the QSS level predicted 
by the Vlasov linear response theory (middle) and the thermal equi- 
librium level with the external field (upper), (a) Short time evolution 
of magnetization, (b) Long time evolution. The horizontal axis is in 
the logarithmic scale, which is scaled as log(Q(r/A') in the inset. 



external field to be small enough. The Vlasov susceptibility is 
computed by taking the limit of h — » 0, and hence the Vlasov 
susceptibility can be computed by the Vlasov linear response 
theory. See Fig. [T]for the high-temperature side. 

The second remark is on the universality of the exponent 
-y^ - p/2. We can estimate D^ for any QSSs by using the 
formula shown in lll3ll . We obtain D^ ~ ^{M)q around 
the critical points for any families of distribution functions /o 
which depend on {q, p) only through the effective one-particle 
Hamiltonian //[/o] and are stable stationary solutions to the 
Vlasov equation. Thus, the exponent ji - f5/2 is obtained for 
such a family of distribution functions. 

The third remark is on the spectrum method to compute sus- 
ceptibility 12211 . The method gives the correct susceptibility in 
the high-temperature side, which is consistent with the present 
result, but does not in the low-temperature side, since the exis- 
tence of integrals is not fully included. The spectrum method 
could be improved by considering the Casimir integrals. 

The fourth remark is that the existence of invariants may 
break some thermodynamic laws. Indeed, local temperature 
in isolated crystalline clusters is not uniform by conservation 
of angular and translational momenta 123 1 . 

We end this article by remarking on the difference between 
dynamics and statistical mechanics. The dynamical system 
gets trapped in QSSs for a long time, therefore, measures on 
experimental setups will show the Vlasov prediction for sys- 
tems with large enough number of particles. 
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with the temperature T - 0.499 < 1/2. The external field 
with h - 0.01 is switched on at f = 0, and the system jumps 
to the QSS predicted by the linear response theory based on 
the Vlasov equation. In the long time regime, the Casimirs 
are no more invariant due to appearance of the collision term, 
and the system goes towards thermal equilibrium. Simula- 
tions indicate that the time scale of relaxation from the QSS 
to thermal equilibrium grows linearly with A^, as found in in- 
homogeneous QSSs without external fields 12 ill . 

Finally five remarks are in order 

The first remark is on the validity of the linear response the- 
ory to discuss critical behaviours. The linear response theory 
assumes that /c - /o is small enough, and it seems not pos- 
sible to discuss divergence of susceptibility at the first sight. 
However, for instance in the low-temperature side, what we 
have to keep is \hx'^\ <s: {M)q. Thus, the linear response the- 
ory may be valid even around the critical point by setting the 
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